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ions,  and/or  findings  contained  in  this  report  are  those  of 
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policy,  or  decision,  unless  so  designated  by  other 
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rain  compression  of  a  rectangular  block  is  used  as  a  model 
vestigate  the  dynamics  of  shear  band  development  from  an 
mogeneity.  The  material  is  characterized  as  a  von  Mises 
plastic  solid,  with  a  hardness  function  that  exhibits  a  local 
ardless  of  whether  the  material  is  hardening  or  softening, 
n  development  involves  the  evolution  of  finger-like  contours 
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Block  20  -  Continuation  of  Abstract 


emanating  from  the  inhomogeneity  at  45  to  the  compression  axis.  Once  a  given 
strain  contour  crosses  the  specimen,  it  fans  out  about  its  initial  direction 
of  propagation.  For  a  softening  solid,  this  fanning  out  ceases  for  some 
strain  level  greater  than  the  strain  at  the  hardness  maximum  and  further 
straining  takes  place  in  an  ever  narrowing  band.  Many  of  the  qualitative 
features  of  shear  band , development  under  dynamic  loading  conditions  are  the 
same  as  under  quasi-static  loading  conditions,  but  a  significant  retardation 
of  shear  band  development  due  to  inertial  effects  is  found. 
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1.  Introduction 


Large  inelastic  deformations  are  frequently  accompanied  by  the  formation  of  intense  bands  of 
localized  shearing.  Often  the  large  localized  strains  in  a  shear  band  precipitate  a  shear  fracture. 
In  other  circumstances  shear  bands  do  not  lead  to  fracture  but  localized  shearing  becomes  an 
important  mechanism  for  subsequent  plastic  deformation.  Hence  shear  bands  have  a  dual 
significance;  as  a  precursor  to  fracture  and  as  a  mechanism  of  large  strain  plastic  response. 

When  the  solid  is  considered  to  be  deforming  quasi-statically  and  isothermally,  shear  band  lo¬ 
calizations  can  be  analyzed  as  a  “material  instability,”  within  a  framework  due  to  Hadamard 
(1903),  Thomas  (1961),  Hill  (1962),  Mandel  (1966)  and  Rice  (1977).  Within  this  framework,  a 
material  element  is  considered  subject  to  all  around  displacement  boundary  conditions  consis¬ 
tent  with  homogeneous  deformations;  different  homogeneous  deformation  states  prevail  outside 
and  inside  the  band  and  the  analysis  determines  the  jumps  in  field  quantities  across  the  band 
interface.  For  rate  independent  solids,  the  onset  of  a  shear  band  localization  can  occur  as  a 
bifurcation  corresponding  to  a  loss  of  ellipticity  of  the  rate  equilibrium  equations,  Hill  (1962), 
Rice  (1977).  For  rate  sensitive  solids  a  bifurcation  is  effectively  excluded  and  the  governing 
rate  equations  remain  elliptic,  although  in  the  presence  of  initial  band  type  imperfections  shear 
band  localizations  do  emerge,  e.g.  Pan  et  al.  (1983). 

The  analysis  of  shear  band  development  in  non-homogeneously  deforming  solids  requires  a 
full  boundary  value  problem  solution.  Such  solutions  have  been  obtained  for  the  quasi-static 
development  of  shear  bands,  see  e.g.  Tvergaard  et  al.  (1981),  Triantafvllidis  et  al.  (1982), 
LeMonds  and  Needieman  (1986)  and  the  review  of  Needieman  and  Tvergaard  (1983),  and  have 
shown  that  shear  bands  do  initiate  with  the  orientation  given  by  the  “material  instability” 
analysis.  Shear  band  propagation,  however,  is  greatly  affected  by  the  inhomogeneity  of  the 
surrounding  fields. 

The  physical  mechanisms  responsible  for  triggering  shear  bands  can  vary  widely.  At  low  strain 
rates,  a  key  feature  of  plastic  material  response  for  localization  is  the  yield  surface  vertex 
structure  implied  by  the  discrete  nature  of  crystallographic  slip.  The  significance  of  a.  vertex 
lies  in  the  reduced  stiffness  to  a  change  in  loading  path  which  permits  shear  bands  to  emerge 
in  strain  hardening  materials,  see  e.g.  Rice  (1977).  At  higher  strain  rates,  thermal  softening 
appears  to  play  a  major  role  in  shear  band  formation,  as  originally  suggested  by  Zener  and 
Holloman  (19-14).  An  overview  of  conditions  for  high  rate  shear  band  formation  in  simple  shear 
and  correlation  with  experiments  is  given  by  Clifton  et  al.  (19X4). 

Profusions  of  internal  shear  bands  are  sometimes  observed,  e.g.  Anand  and  Spitzig  ( 19X0),  and 
these  bands  appear  to  originate  at  internal  inhomogeneities.  Analyses  of  shear  localization  from 
internal  strain  concentrators  have  been  carried  out  by  Abeyaratne  and  Triantafvllidis  (19X1), 
1  vergaard  (19X2,  19X7),  Freund.  Wu  and  Toulins  (1!)X.">)  and  LeMonds  and  Needieman  (19X6). 


Abeyaratne  and  Triantafyllidis  (1981)  investigated  shear  band  initiation  in  quasi-statically  de¬ 
forming  hypo-elastic  and  hyper-elastic  solids  in  plane  strain.  LcMonds  and  Needleman  (198G) 
also  considered  plane  strain  quasi-static  deformations  but  analyzed  the  course  of  shear  band 
development  in  elastic- viscoplastic  solids,  accounting  for  strain  hardening,  rate  dependence, 
thermal  softening  and  heat  conduction.  Tvergaard  (1982,  1987)  has  studied  the  quasi-static 
response  of  porous  plastic  solids  having  two  size  scales  of  void  nucleating  particles  and  where 
early  nucleation  at  relatively  large  particles  triggers  localization  into  a  void  sheet.  In  Freund, 
Wu  and  Toulios  (1985)  the  dynamics  of  shear  band  propagation  was  analyzed  for  a  nonlinear 
elastic  solid  in  anti-plane  strain. 

Here,  t\_>  dynamics  of  shear  band  development  from  internal  material  inhomogeneities  under 
plan  train  compressive  loading  is  used  as  a  model  problem  for  investigating  shear  band  prop¬ 
agation  in  a  rate  sensitive  solid  under  plane  strain  conditions.  The  material  is  characterized  as 
an  elastic-viscoplastic  von  Mises  solid.  In  order  to  simulate  the  phenomenon  of  thermal  soften¬ 
ing  due  to  adiabatic  heating,  the  material  is  taken  to  be  strain  softening.  Although  attention 
is  focussed  on  circumstances  where  such  softening  is  responsible  for  shear  band  development,  it 
is  important  to  emphasize  that  for  more  general  constitutive  descriptions  of  plastic  flow,  strain 
softening  is  not  necessary  to  precipitate  a  shear  band  localization,  Rudnicki  and  Rice  (1975), 
Rice  (1977).  Due  to  the  viscoplastic  formulation  wave  speeds  remain  real  even  in  the  softening 
regime;  in  fact,  regardless  of  the  plastic  pre-strain,  wave  speeds  are  set  by  the  elastic  moduli. 
Furthermore,  as  a  consequence,  the  pathological  mesh  size  effects  encountered  in  localization 
analyses  for  rate  independent  solids  are  eliminated,  Needleman  (1987). 

2.  Problem  Formulation 

The  analysis  is  based  on  a  convccted  coordinate  Lagrangian  formulation  of  the  field  equations 
with  the  initial  unstressed  state  taken  as  reference.  All  field  quantities  are  considered  to  be 
functions  of  convccted  coordinates,  y',  which  serve  as  particle  labels,  and  time  t.  The  position, 
relative  to  a  fixed  Cartesian  frame,  of  a  material  point  in  the  initial  configuration  is  denoted 
by  x.  In  the  current  configuration  the  material  point  initially  at  x  is  at  x.  The  displacement 
vector  u  and  the  deformation  gradient  F  are  defined  by 

Ox 

u  —  x  —  x  F  =  —  (2.1) 

Ox 

Base  vectors  in  the  reference  configuration  (unbarred)  and  in  the  current  configuration  (barred) 
are  given  by 


where  g,J  and  glJ  are,  respectively,  the  inverses  of  the  metric  tensors  ytJ  -  gj  -gj  and  Oij  —  g,  gj. 


The  momentum  balance  can  be  expressed  either  in  terms  of  the  symmetric  Cauchy  stress  tensor 
a  or  the  nonsymmetric  nominal  stress  tensor  s.  These  are  related  to  the  force,  df,  transmitted 
across  a  material  element  by 


df  =  v -<rdS  =  v  •  sdS  (2.4) 

Here,  dS  and  u  give  the  area  and  orientation,  in  the  current  configuration,  of  a  material  element 
that  had  area  dS  and  orientation  v  in  the  reference  configuration  and  the  stress  measures  s  and 
a  are  related  by 


s  =  F  1  r  t  -  det(F)a 


The  dynamic  principle  of  virtual  work  is  written  as 


/  sji6Fi:idV=  [rtuidS-  /  p^- 6 utdV  (2.6) 

JV  Js  Jv  vl- 

where  V ,  S  and  p  are  the  volume,  surface  and  mass  density,  respectively,  of  the  body  in  the 
reference  configuration  and 


T  =  s3iUj  (2.7) 

The  specific  problem  considered  is  plane  strain  compression  of  an  initially  square  block  of 
dimension  ‘2Iiq  X  2/iq.  The  i/1  —  y2  plane  is  taken  to  be  the  plane  of  deformation  and  the 
compression  axis  is  aligned  with  the  y 2  direction.  The  origin  of  the  coordinate  system  is  placed 
at  the  center  of  the  block  and  identical  displacement  histories  are  imposed  at  y 2  =  ±/i0  and 
the  sides  yl  =  ±/io  remain  traction  free.  Attention  is  restricted  to  deformations  that  remain 
symmetric  about  both  i/1  =  0  and  y 2  —  0.  The  boundary  conditions  for  the  quadrant  analyzed 
numerically  are 


7/ 1  =  0  /'•  =  0  at  i/  =  0 
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where  (')  is  d(  )/dt. 

The  function  U(t)  in  (2.10)  is  taken  as 


U(t)  = 


jVrit/h), 

l  Vi, 


for  t  <  1 1 ; 
for  t  >  t\ 


(2.12) 


with  Vi  being  a  prescribed  velocity. 

The  material  is  characterized  as  an  elastic-viscoplastic  von  Mises  solid.  The  total  rate  of 
deformation,  D,  is  written  as  the  sum  of  an  elastic  part,  De,  and  a  plastic  part  Dp,  with 


De 


1  +  //.  v 

E  T  ~  £  r  "  11 


(2.13) 


=  =  (2.14) 

where  f  is  the  Jaumann  rate  of  Kirchholf  stress,  I  is  the  identity  tensor,  f  :  I  is  the  trace  of  t, 
1  is  the  effective  plastic  strain  rate,  E  is  Young’s  modulus,  v  is  Poisson’s  ratio  and 

T'  =  t-^(t:I)I  (2.15) 

Combining  (2.13)  and  (2.14)  and  inverting  to  express  f  in  terms  of  D  gives  a  relation  of  the 
form 


r=L:D-P  P  =  <L:p  =  cp:L 


(2.  Hi) 


where  L  is  t lie*  tensor  of  elastic  moduli  and  the  dyadic  product  L  :  p  denotes  pik. 


On  the  current  base  vectors  (2.1G)  has  the  form 

Tij  =  Liikt Etk  -  P”  (2.17) 

The  Lagrangian  strain  rate  appears  in  (2.17)  via  the  identity  Ejj  =  g,  Dgj  and  its  components 
on  the  reference  base  vectors  are  given  by 

Eij  =  ^  (FkFkj  +  FkFk)  (2.18) 

The  constitutive  relation  (2.17)  can  be  expressed  in  terms  of  the  convected  derivative  of  Kirch- 
hoff  stress  on  the  current  base  vectors,  i.e.  r,J  =  g '  t  ■  gJ,  as 

fij  =Cijk,Elk-  Pij  (2.19) 

where 

Cijk‘  =  L'-kl  -  i  gikrjl  +  (jurjk  +  gjkTil  +  c/rik  (2.20) 

The  uniaxial  response  is  represented  by  the  power  law  rate  relation 

*  =  *o(*/fl(0]1/w  (2.21) 

The  hardness  function  g((),  with  £  =  /  f (U,  represents  the  effective  stress  versus  effective  strain 
response  in  a  tensile  test  carried  out  at  a  strain-rate  such  that  (  =  co-  In  (2.21 ).  an  is  a  reference 
strength  and  m  >  0  is  the  strain  rate  hardening  exponent. 

On  quite  general  grounds,  for  metals  undergoing  plastic  deformation  due  1o  dislocation  motion, 
the  plastic  strain  rate  is  appropriately  taken  as  a  function  of  the  stress  and  the  current  state. 
Rice  ( 1970).  The  relation  (2.21)  employs  the  idealization  of  a  single  parameter  characterization 
of  the  current  state.  A  fundamental  distinction  between  rate  dependent  material  behavior  as 
embodied  in  (2.21 )  and  rate  independent  behavior  is  that  the  plastic  strain  rate  in  (2.21 )  does 
not  depend  on  incremental  quantities.  For  structural  metals  at.  room  temperature,  value's  of  in 
in  the  range  0.002  to  0.02  are  representative;  see.  for  example,  data  collected  by  (Ihosh  ( 1977). 
liven  though  there  is  no  explicit  plastic  yielding  condition  in  the  constitutive  formulation,  for 
this  range  of  m,  there  is  an  apparent  rate  dependent  yield  strength  and  an  “unloading"  like 
response  on  strain  rate  reversal,  see  e.g.  Needleman  (19X7). 

■  ) 


The  power  law  description  of  material  strain  rate  sensitivity  in  (2.21)  is  not  a  fundamental 
physical  relation,  but  it  can  provide  a  reasonable  representation  of  the  strain  rate  sensitivity 
of  common  structural  metals  over  several  decades  of  strain  rate.  The  physical  background  for 
power  law  viscoplasticity,  and  its  limitations,  are  discussed  by  Asaro  (1983)  and  Klopp,  Clifton 
and  Shawki  (1985). 

The  hardness  function  g(l)  is  taken  to  have  the  form 


'j(i)  =  tfo 


(1  +  c/f0)v 
1  +  (f/fi)2 


(2.22) 


For  —  oo,  (2.22)  reduces  to  a  power  law  hardening  relation.  As  ci  decreases,  g(i)  exhibits 
strain  softening  at  decreasing  values  of  c.  Figure  1  shows  curves  of  g(c)  versus  c  with  N  =  0.1 
and  c0  =  0.00218  and  for  ex  =  104  x  c0,  200  x  c0  and  100  x  c0.  The  finite  element  calculations 
to  be  discussed  subsequently  are  based  on  the  three  hardness  functions  shown  in  Fig.  1. 

Since  the  plastic  strain  rate  in  (2.21)  depends  on  the  current  state  and  not  on  incremental 
quantities,  the  propagation  of  disturbances  depends  on  the  elastic  stiffnesses,  which  is  consis¬ 
tent  with  observation,  see  e.g.  Clifton  (1972).  In  particular,  the  speed  at  which  incremental 
disturbances  propagate  does  not  depend  explicitly  on  whether  the  material  is  hardening  or 
softening.  Over  the  range  of  stress  and  deformation  states  arising  in  the  numerical  examples 
here,  the  maximum  ||  r  ||<  O.Ul  x  E  and  maximum  principal  logarithmic  strains  less  than  0.5, 
wave  speeds  remain  nearly  constant  at  their  linear  elastic  values,  even  in  the  softening  regime. 

3.  Numerical  Method 

The  discretization  is  carried  out  in  the  reference  configuration  with  the  convected  coordinates, 
y * ,  being  the  independent  variables.  Linear  displacement  triangular  elements  arc  used  arranged 
in  a  quadrilateral  pattern  consisting  of  ‘'crossed"  triangles.  The  "crossed"  triangular  quadrilat¬ 
erals  are  used  both  because  of  their  ability  to  represent  incompressible  deformations  without 
mesh  locking,  Nagtcgaal  et  al.  ( 1971),  and  because  of  their  ability  to  reproduce  localized  defor¬ 
mation  modes  that  are  oblique  to  the  quadrilateral  mesh.  Tvergaard  et  al.  (19X1),  Needleman 
and  Tvergaard  (1983).  A  lumped  mass  matrix  is  used  instead  of  the  consistent  mass  matrix 
in  (3.6),  since  this  has  been  found  preferable  for  explicit  time  integration  procedures,  from  the 
point  of  view  of  accuracy  as  well  as  computational  efficiency,  Krieg  and  Key  ( 1973).  When  the 
finite  element  approximations  of  the  displacement  components  u‘  are  substituted  into  (2.6). 
and  the  integrations  are  carried  out.  the  discretized  equations  ol  motion  are  obtained 


where  M  is  a  mass  matrix,  U  is  a  nodal  displacement  vector  and  R  is  a  nodal  force  vector. 


An  initial  inhomogeneity  is  specified  that  corresponds  to  a  local  reduction  in  flow  strength, 
in  (2.1 1 )  of  the  form 


°o  (<j  <y')  =  °o 


1  -  oexp{[(i/‘  -  xjI)2  +  (y2  -  J/o )']/ro } 


(3.2) 


where  j/q  and  Uo  specify  the  location  of  the  inhomogeneity  and  ro  and  a  give  the  inhomogeneity 
size  and  magnitude,  respectively. 

At  t  =  0.  the  rectangular  block  under  consideration  is  taken  to  be  stress  free  and  motionless.  The 
equations  of  motion  (3.1)  are  integrated  using  a  central  difference  scheme  that  can  be  regarded 
as  a  member  of  the  Newmark  family  of  time  stepping  algorithms,  Hughes  and  Belytschko  ( 1983). 
At  each  time  increment  the  stress  state  is  updated  using  the  rate  tangent  method  of  Peirce  et 
al.  (1984).  This  combination  of  a  central  difference  time  stepping  scheme  in  conjunction  with 
the  tangent  modulus  stress  updating  algorithm  is  discussed  by  Moran  (1987)  and  has  also  been 
used  in  Needleman  (1987)  and  Tvergaard  and  Needleman  (1987).  For  the  present  problem,  the 
basic  equations  take  the  form 


Un+1  =  Un  +  A<Vn+i(A/)2^ 

2  at 


(3.3) 


flVn+l 

Ot 


=  -M-,R„+i 


(3,1) 


V„+, 


OVn 

Ot 


+ 


Wn  + 1 
Ot 


(3.5) 


Here,  V  is  the  vector  of  nodal  displacement  and  velocity,  M  1  denotes  the  inverse  of  the  mass 
matrix  and  the  subscripts  n  and  n  +  1  refer  to  quantities  evaluated  at  /„  and  /,1  +  ] .  respectively. 

At  each  integration  point,  the  stress  tensor  is  updated  via 


'  N  +  l 


+  A/f‘ 


( 3  .(> ) 


1  li,  stress  rate  components.  f,J ,  are  computed  from  (2.19)  and  (2.20)  with  L  and  P  replaced 
by  the  corresponding  rale  tangent  quantities. 


L  =L-rr^(L:p)Q(L:p)  p  =tt^L:p 


where  A  ®  B  denotes  the  tensor  product  A'i Bkl  and 


ft  =  eo[d/y(c)]1/m  i-{6At)h-^ 


( on  ( dt 

,  =  p;L:p-  *  * 


In  the  second  of  (3.8),  0  is  a  parameter  of  the  integration  scheme  that  is  taken  to  be  unity  in 
the  numerical  calculations  here. 

The  effective  plastic  strain  is  updated  using  +  ]  =  ln  +  Ati  with 


l  =  —  =  + - - - (L  :  p)  :  D 

At  1+^(1  +  0* 


(3.10) 


The  constitutive  updating  calculations  use  (Un+i  -  U n)/At  to  represent  the  displacement  rate 
components  in  (2.18). 

In  the  limit  £  —  oc,  (3.7)  through  (3.10)  reduce  to  their  rate  independent  counterparts.  On 
the  other  hand,  for  £  =  0,  which  results  from  setting  0  =  0  in  the  second  of  (3.8),  the  Euler 
integration  scheme  is  recovered. 

4.  Results 

In  the  calculations  here,  all  material  parameters  are  kept  fixed  except  for  fi  which  governs  the 
rate  of  softening.  Although  the  aim  here  is  not  to  model  any  particular  material,  the  fixed 
property  values  are  taken  to  be  representative  of  structural  steels;  E  =  211  (5 Pa.  rr0  =  1(>0 
MPa.  (0  =  d0/£  =  0.00218,  u  =  0.3.  .V  =  0.1,  m  =  0.01  and  c0  =  0.002  sec-1.  The 
initial  inhomogeneity  in  (3.2)  is  specified  by  y$  =  1/5  =  0,  a  =  0.2  and  i'o  =  0.1  mm.  With 
p  -  7800  kg/m3  =  7.8  X  10-3MPa/(m/sec)‘,  the  elastic  wave  speeds  are 


r»  =  y  —  =  3.23  X  10  in /sec 


1 2(7 ( 1  -r)/(  1  -  2/0 


=  (>.03  x  lO'm/sec  (1.1) 


where  2( !  =  /./( 1  +  //). 


Three  values  of  the  imposed  velocity  are  considered;  V'i  --  -3  m/sec,  Vj  =  -10  m/sec  and 
Vj  =  -30  m/sec,  which  correspond  to  |  Vj  |  /co  =  0.920  x  10-3,  |  \\  |  /c0  =  3.10  X  10-3  and 
I  Vj  |  /c0  =  9.29  X  1 0— 3 ,  respectively.  In  each  case  the  rise  time,  in  (2.12),  is  0.1  X  10_f> 
sec.  Also,  the  block  size  ho  is  taken  to  be  1  mm,  so  that  the  nominal  strain  rates  for  the  three 
imposed  velocities  are  3  X  103,  10*  and  3  X  1  O'* ,  respectively. 

A  uniform  finite  element  mesh  is  used,  with  the  elements  arranged  so  that  at  a  strain  a  little 
larger  than  when  the  peak  hardness  is  attained,  the  element  diagonals  are  oriented  for  45°  shear 
bands.  For  the  majority  of  calculations,  with  ci  =  200  X  fo,  a  26  X  20  element  mesh  is  used  (the 
deformed  mesh  is  shown  subsequently  in  Fig.  8),  while  with  C]  =  100  x  c0  a  24  x  20  mesh  is 
used.  In  each  case  the  20  element  division  is  in  the  y~  direction.  A  variable  time  step  method 
is  used  with  one  criterion  limiting  the  time  step  being  the  requirement  that  the  time  step  be  no 
greater  than  the  time  for  a  wave  travelling  at  velocity  C[  to  traverse  half  the  smallest  element 
side.  With  the  ratio  of  wave  speeds  in  (4.1),  this  is  about.  94%  of  the  time  for  a  wave  travelling 
at  the  shear  wave  speed  c0  to  traverse  the  smallest  element  side.  Other  time  step  limits  are  set 
by  stability  requirements  for  the  rate  tangent  modulus  integration  of  the  constitutive  relations. 
However,  in  the  problems  considered  here  these  are  rarely  operative  and  the  time  step  limits 
set  by  stability  requirements  for  the  explicit  integration  of  the  equations  of  motion  determine 
the  maximum  permissible  time  step. 

Figure  2  shows  curves  of  average  stress  versus  end  displacement,  with  an  imposed  velocity  of  10 
m/sec,  for  each  of  the  hardening  relations  in  Fig.  1.  The  average  stress  (per  unit  out-of-plane 
thickness)  is  given  by 


where  li  is  the  current  width  at  y2  =  /i0  and  the  end  displacement  is  U\  =  \\<lt. 

The  initial  response  is  identical  for  the  three  hardening  relations.  Once  extensive  plastic  defor¬ 
mation  occurs,  dynamic  ellect.s  are  rapidly  damped.  For  the  softening  solids  there  is  an  abrupt 
stress  drop  that  is  associated  with  shear  band  formation. 

Before  considering  - 1 1« -,i i  band  piop.ng.it  am.  tin*  dynamics  ol  strain  development  for  the  hard¬ 
ening  solid  F  .'hown  in  fig.  3.  Sti.iin  fonioiiis  oiiginate  at  the  inhomogeneit  v  and  evolve  as 
linger -like  ronton  t  >  at  I  V  to  t  In-  *  m; ;  ■;  •  --.oil  axis.  Ome  a  given  contour  crosses  the  specimen, 
it  Ian-  out  a  bou  t  its  initial  * ;  i :  >  .  j  u  op.igal  ion .  1  n  1  ig.  3.  the  pro  pa  gat  ion  speed  along 

the  br  directum  is  v . > t )  m  -.■<  jot  t|  ••  o.os  ,,.n'our  and  1  lilt)  m/sec  lor  the  0.09  contour. 

I  or  t  In*  soil  ening  -olid  i  a  I  ig.  I .  - :  i . . 1  •  \  •  I  <  >  p  1 1  n  - 1 1 1  mil  i.  ill  v  l<  4  lows  the  pattern  shown  in  Fig.  3. 

A  shear  baud  biiut  <  at  um  o.  i  1 1 1  -  t||i*  u  tide  i  Is  i  tig  t  ate  independent  solid  at  t  lie  at  t  a  in  men!  of 


the  maximum  hardness  and  the  45°  shear  band  orientation  is  the  most  critical  one.  Eventually, 
for  a  strain  level  beyond  the  peak  in  Fig.  1,  spreading  in  the  direction  perpendicular  to  the 
initial  propagation  direction  slows  down  and  further  straining  occurs  in  a  band.  The  ultimate 
course  of  events  is  essentially  determined  at  the  stage  of  loading  shown  in  Fig.  4b;  continually 
increasing  strain  in  a  shear  band  with  very  little  change  in  strain  outside  the  band.  Between 
Figs.  4a  and  4b  the  speeds  of  propagation  of  the  0.07  and  0.0S  contours  are  2200  m/sec  and 
1000  m/sec,  respectively.  Between  Figs.  4b  and  4c,  the  propagation  speed  of  the  0.09  contour 
along  the  band  is  2900  m/sec,  which  is  approaching  the  shear  wave  speed.  For  the  case  in  Fig. 
2  with  Ci  =  200  X  to,  plastic  strain  contour  speeds  at  a  similar  stage  of  band  development  are 
somewhat  slower,  around  1300  m/sec. 

Figure  5  shows  the  effect  of  prescribed  velocity  oil  the  average  stress  versus  end  displacement 
response  for  a  material  characterized  by  the  hardness  function  with  tj  =  200  X  (q.  With 
V’l  =  —3  m/sec,  corresponding  to  an  average  strain  rate  of  3000/sec.  the  response  is  essentially 
quasi-static  as  soon  as  there  is  any  significant  plasticity.  When  the  velocity  is  increased  to 
-10  m/sec,  (this  is  the  same  case  as  B  in  Fig.  2)  dynamic  effects  persist  much  longer,  but 
have  been  essentially  damped  out  by  the  time  localization  initiates.  At  \\  —  -30  m/sec, 
significant  inertial  effects  are  evident  through  localization  initiation.  If  the  imposed  velocity 
were  increased  much  beyond  -30  m/scc,  the  average  stress  would  change  sign  at  the  first 
trough,  indicating  a  separation  event,  subsequent  to  which  imposition  of  (2.11)  would  not  be 
the  physically  appropriate  boundary  condition. 

Quasi-static  calculations  were  carried  out  using  the  same  material  properties  and  using  an 
identical  finite  element  mesh  for  the  cases  in  Fig.  5  with  \'i  =  -3  m/sec  and  li  =  -30 
m/scc.  With  V\  =  -3  m/scc.  except  for  the  lack  of  “wiggles."  the  quasi-statically  computed 
stress  strain  curve  coincides  with  that  in  Fig.  5  up  to  the  drop  due  to  shear  localization. 
The  quasi-static  calculation  gives  localization  at  a  strain  0.0000  smaller  than  given  by  the 
dynamic  calculation  in  Fig.  5,  but  other  than  this  small  shift,  the  quasi-static  post-localization 
response  is  indistinguishable  from  the  dynamic  response  shown.  The  case  analyzed  here  is 
one  for  which  there  are  equilibrium  solutions  with  increasing  end  displacement  during  shear 
band  development.  In  some  circumstances,  see  e.g.  Tvergaard  et  al.  (19X1),  there  may  be  no 
equilibrium  solutions  having  increasing  end  displacement  as  localization  proceeds  and  in  such 
cases  dynamic  effects  may  be  significant  for  localization  development  even  at  very  low  imposed 
velocities. 

In  the  quasi-static  analysis  with  a  prescribed  velocity  of  —39  m/sec.  the  order  of  magnitude 
increase  in  the  nominal  applied  strain  rate,  from  3  X  It)’  to  3  X  It)1,  gives  an  increase  in  stress 
level  of  about  '1%  (  If)0  01  as  1.02),  but  shear  localization,  as  signalled  by  the  drop  in  the  overall 
stress-strain  curve,  occurs  at  nearly  the  same  end  displacement  value  as  with  1  |  —  —3  m/sec. 
1  his  behavior  is  expected,  since  analvses  of  the  effect  of  rate  sensitivity  on  plastic  instabilities 


10 


in  rigid-plastic  power  law  rate  hardening  materials,  e.g.  Hutchinson  and  Neale  (1977),  show 
that  while  the  instability  strain  depends  sensitively  on  the  strain  rate  hardening  exponent,  it 
is  independent  of  the  rate  Hence,  the  delay  in  shear  band  initiation  with  increasing  imposed 
velocity  evident  in  Fig.  5  is  a  consequence  of  material  inertia  and  not  of  material  strain  rate 
sensitivity. 

There  is  a  significant  distinction  between  the  dynamic  and  quasi-static  formulations  regarding 
characteristic  time  scales.  In  the  quasi-static  formulation  there  is  a  single  material  time  scale, 
1/eo,  while  in  the  dynamic  formulation  there  are  two;  1/co  and  a  problem  dependent  time  scale 
that  is  a  wave  speed  divided  by  a  characteristic  length,  e.g.  co//io-  The  rate  of  strain  develop¬ 
ment  scales  with  Vi/(oho  in  the  quasi-static  formulation,  but  not  in  the  dynamic  analyses,  due 
to  inertial  effects. 

Figures  6  and  7  show  the  course  of  shear  band  development  for  V\  =  —3  m/sec  and  Vj  =  —30 
m/sec,  respectively.  There  is  a  clear  effect  of  imposed  velocity  on  the  speed  of  propagation 
of  the  strain  contours.  For  example,  in  Fig.  6  the  0.15  contour  propagates  along  the  band 
at  a  speed  between  590  m/sec  and  080  m/sec,  while  in  Fig.  7  the  corresponding  propagation 
speed  is  2500  m/sec.  The  qualitative  features  of  shear  band  development  do  not  depend  on  the 
imposed  rate,  but  it  is  evident  in  Figs.  6  and  7,  as  well  as  in  Fig.  5,  that  there  is  a  delay  in 
shear  band  initiation  with  increasing  rate. 

Figure  8  shows  deformed  finite  element  meshes  at  two  stages  of  deformation.  Fig.  8a  corre¬ 
sponds  to  the  stage  of  deformation  in  Fig.  7c.  Even  though  the  shear  band  has  already  formed, 
this  is  difficult  to  ascertain  from  the  deformed  mesh.  At  the  stage  of  deformation  in  Fig.  8b,  the 
mesh  is  incapable  of  resolving  the  sharp  strain  gradients  across  the  band  and  the  discreteness 
of  the  mesh  significantly  affects  the  computed  response.  This  figure  is  included  to  illustrate  the 
ability  of  the  mesh  to  develop  a  narrow  shear  band,  which  confirms,  for  example,  that  the  initial 
thickness  of  the  band  evident  in  Figs.  3,  -1,  G  and  7  is  an  actual  reflection  of  the  inhomogeneity 
size  and  is  not  a  mesh  artifact. 

In  order  to  explore  the  effect  of  mesh  on  the  computed  response,  two  additional  meshes  were 
used  for  the  case  with  V\  =  — 10  m/sec  and  fi  =  100  X  c0-  !>'  one  ease  the  mesh  aspect  ratio  was 
varied,  while  in  the  other  case  the  aspect  ratio  was  kept  fixed,  but  a  finer  mesh  was  used.  The 
computed  average  stress  versus  end-displacement  curves  are  shown  in  Fig.  9.  In  the  29  x  20 
mesh  the  band  is  spread  across  two  rows  of  quadrilaterals.  'This  gives  rise  to  a  significant  mesh 
induced  delay  in  the  onset  of  shear  band  formation  even  though  this  mesh  has  more  elements 
than  the  reference  mesh.  Increasing  the  number  of  elements,  but  keeping  a  fixed  aspect  ratio 
leads  t  o  a  slight  ly  earlier  onset  of  local  izat  ion.  which  is  in  part  due  to  a  bet  ter  resolution  of  the 
strain  concent  ration  at  the  inhomogeneity.  The  liner  mesh  can  also  resolve  larger  gradients  and 
so  exhibits  a  greater  stress  drop  before  mesh  induced  stillening  occurs,  flic  post -localization 


rise  in  the  stress  versus  end-displacement  curve  for  the  29  X  20  mesh  is  a  consequence  of  mesh 
induced  stiffening  when  gradients  in  the  band  become  very  large.  A  similar  stiffening  occurs 
for  the  other  two  meshes  in  Fig.  9,  but  occurs  later  in  the  deformation  history  with  improved 
resolution  of  the  band. 

The  speed  of  propagation  of  the  constant  plastic  strain  contours  does  not  differ  much  between 
these  three  meshes.  For  the  24  X  20  mesh  in  Fig.  4,  the  propagation  speed  along  the  band  of 
the  0.08  plastic  strain  contour  is  1600  m/sec.  At  a  similar  stage  of  deformation  with  the  36  X  30 
mesh  the  0.08  contour  is  propagating  at  1900  m/sec  and  at  1600  m/sec  with  29  x  20  mesh.  The 
difference  in  propagation  speed  found  for  the  36  x  30  mesh  is  attributable  more  to  variations 
in  speed  as  the  contour  propagates  along  the  band  than  to  the  increased  resolution;  with  the 
24  x  20  mesh  in  Fig.  4,  a  propagation  speed  of  2200  m/sec  is  found  for  the  0.07  contour. 

5.  Discussion 

Key  features  of  the  phenomenology  of  shear  band  development  under  dynamic  loading  condi¬ 
tions  are  the  same  as  under  quasi-static  loading  conditions.  Shear  bands  develop  at  45°  from 
the  compression  axis,  as  expected  based  on  a  quasi-static  “material  instability”  analysis.  Under 
quasi-static  loading  conditions,  for  the  viscoplastic  solids  considered  here,  the  strain  at  which 
shear  bands  initiate  is  an  outcome  of  the  competition  between  the  destabilizing  influence  of 
the  inhomogeneity  and  the  stabilizing  efTect  of  material  strain  rate  sensitivity.  Under  dynamic 
loading  conditions,  a  delay  in  shear  band  development  arises  due  to  inertial  effects.  Also,  since 
dynamic  effects  induce  non-homogeneous  deformations,  an  inhomogeneity  is  not  required  to  ini¬ 
tiate  localization  but,  in  the  circumstances  analyzed  here,  the  inhomogeneity  plays  the  major 
role  in  triggering  the  onset  of  shear  localization. 

The  key  to  the  development  of  shear  bands  in  the  calculations  here  is  the  speed  of  propagation 
of  plastic  strain  contours  in  the  direction  normal  to  the  eventual  band.  This  is  consistent 
with  what  would  be  expected  based  on  a  rate  independent  shear  band  bifurcation  analysis, 
where,  at  bifurcation,  the  propagation  velocity  of  shear  waves  normal  to  the  band  vanishes. 
The  vanishing  of  this  wave  speed  is  also  the  shear  band  initiation  criterion  that  emerges  from 
the  perturbation  analysis  of  Freund,  \Yu  and  Toulios  (1985).  The  process  of  band  formation,  a 
cascade  of  fmger-like  plastic  strain  contours  emanating  from  the  inhomogeneity  and  from  a  small 
region  where  the  band  intersects  the  impacted  surface,  in  a  continually  sharpening  band,  is  one 
that  makes  defining  a  band  propagation  speed  somewhat  arbitrary.  Although  there  is  no  single 
event  that  can  be  identified  with  shear  band  initiation,  there  is  a  reasonably  well  identifiable 
strain  in  Figs.  4,  6  and  7  above  which  all  further  straining  occurs  within  a  narrowing  band.  A 
band  propagation  speed  can  be  identified  with  the  propagation  speed  ol  some  particular  strain 
contour  along  t lie  band,  but  this  speed  is  not  associated  with  some  sharp  gradient  propagating 
through  the  material. 
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The  main  strain  gradients  occur  across  the  band  and  large  strain  gradients  do  not  develop  along 
most  of  the  band;  substantial  strain  gradients  in  the  band  direction  are  confined  to  very  near 
the  inhomogeneity  and  to  a  small  region  where  the  band  intersects  the  impacted  surface.  This 
can  be  seen  in  Figs.  4,  6,  and  7;  for  example,  even  in  Fig.  7c,  the  element  in  the  center  of 
the  inhomogeneity  has  an  effective  plastic  strain  of  0.33.  while  in  the  adjacent  element  in  the 
shear  band  the  effective  plastic  strain  value  is  0.25  (clearly,  the  mesh  used  here  does  not  resolve 
the  details  of  field  distributions  around  the  in  homogeneity).  The  continually  steepening  strain 
gradient  across  the  band  is  much  like  that  in  one  dimensional  simple  shear  calculations,  e.g. 
Wada  et  al.  (1978),  Wu  and  Freund  (1984),  Shawki,  Clifton  and  Majda  (1987)  and  Needleman 
(1987).  Such  narrowing  shear  bands  are  observed,  Marchand  and  Duffy  ( 1987).  In  actuality,  of 
course,  under  the  conditions  of  temperature,  strain  and  strain  rate  that  prevail  in  a  narrow  shear 
band,  material  behavior  is  not  described  by  the  sort  of  simple  phenomenological  constitutive 
relation  used  in  this  investigation.  Also,  some  material  dependent  length  scale  ultimately  sets 
a  minimum  band  thickness. 

The  situation  analyzed  in  this  investigation  is  one  where,  just  prior  to  shear  band  formation, 
the  stress  and  deformation  fields  are  rather  uniform.  When  shear  band  development  takes  place 
from  a  more  non-uniform  deformation  state,  the  gradients  in  the  background  field  retard  shear 
band  growth,  as  is  shown  in  quasi-static  analyses,  Tvergaard  et  al.  (1981)  and  Trianta fyllidis 
et  al.  (1982).  Under  dynamic  loading  conditions,  when  there  are  strong  unfavorable  gradients, 
shear  bands  are  likely  to  propagate  in  a  more  crack-like  fashion  than  in  the  cases  analyzed  here. 

A  strain  softening  constitutive  relation  has  been  used  as  a  simple  model  for  a  thermal  softening 
solid.  For  an  actual  thermal  softening  solid,  since  the  softening  behavior  is  governed  by  the 
current  rate  of  plastic  dissipation,  increased  softening  would  accompany  an  increase  in  imposed 
strain  rate  and  this  would  tend  to  counteract  the  delay  in  localization  due  to  inertial  effects. 
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Figure  Captions 

Figure  1  Curves  of  hardness  y  versus  effective  plastic  strain,  t.  from  (2.22)  with  =  0.1, 
c0  =  0.00218  and  (A)  c,  =  10*  X  c0,  (U)  >x  =  200  x  c0  and  ('/)  c,  =  100  x  r0. 

Figure  2  Average  stress  versus  end  displacement  curves  with  a  prescribed  velocity  of  l'i  =  —10 
m/sec.  The  fixed  material  properties  and  initial  inhomogeneity  are  specified  in  t lie 
text  and  curves  (A),  (11)  and  (C)  correspond  to  the  three  hardness  functions  shown  in 
Fig.  1. 

Figured  Contours  of  constant  effective  plastic  strain,  c,  for  a  hardening  solid.  ct  =  10' xc0.  The 
fixed  material  properties  and  initial  inhomogeneity  are  specified  in  the  text  and  the 
prescribed  end  velocity  is  l't  =  — 10  m/sec;  (a)  t  =  0.14  x  10_(’  sec  (l>)  t  =  0.71  x  10_<i 
sec  (c)  I  =  0.9S  x  10_l>  sec  (d)  t  =  7.25  X  10~°  sec. 

Figure  d  Contours  of  constant  effective  plastic  strain,  c,  for  a  softening  solid,  t  ]  =  100  xc0.  The 
fixed  material  properties  and  initial  inhoinogeneitv  are  specified  in  the  text  and  the 
prescribed  end  velocity  is  l’i  =  -10  m/sec;  (a)  t  =  1.75  X  10_G  sec  (b)  t  =  5.02  X  10-G 
sec  (c)  t  =  5.28  x  L0~6  sec. 

Figure  5  Average  stress  versus  end  displacement  curves  for  a  softening  solid.  <,  =  200  x  r0  with 
(A)  I)  -  —  00  m/sec  (11)  l'j  -  -10  m/sec  and  (C)  \’i  =  -d  m/sec. 

Figure  0  Contours  of  constant  effective  plastic  strain.  < .  for  a  softening  solid,  <  |  =  200  X  <0.  The 
fixed  material  properties  and  initial  inhomogeneity  are  specified  in  the  text  and  the 
prescribed  end  velocity  is  l'i  =  -d  m/sec;  (a)  t  -  28.07 x  10-1’  sec  (b)  /  =  20. H) x  K)'1' 
sec  (c)  I  =  29.71  X  10-'1  sec. 

Figure  7  ( 'ontours  of  constant  effective  plastic  strain,  r.  for  a  softening  solid,  < ,  =  200  x  <()-  The 
fixed  material  properties  and  initial  inhoinogeneitv  are  specified  in  the  text  and  the 
prescribed  end  velocity  is  Fj  =  -.00  m/sec;  (a)  t  =  d.2  1  x  Id-1’  sec  (!>)/  =  0.50  x  111-'' 
sec  (c)  /  =  .0.70  X  10  '*  s(<c. 

Figure  X  Deformed  finite  element  meshes  for  a  softening  solid.  (  i  -  200  Xm.  F.ach  (|iiadrilat<‘ral 
consists  of  four  "crossed"  triangular  elements.  The  fixed  material  properties  and  initial 
inhomogeneity  are  specified  in  the  text  and  the  prescribed  end  velocity  is  1)  -  -00 
m/sec;  (a)  /  -  0.70  x  10_,>  sec  (!>)/-  !.7!)  x  It)-1'  sec. 

figure  !)  Stress  versus  end  displacement  for  a  softening  solid.  (|  -  200  x  with  a  prescribed 
velocity  l  i  -  10  m/sec  using  t  liree  meshes;  j  ,\  )  g"l  x  20  (  111  2  I  x  20  and  ( ( ' )  00  x  00. 


